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Summary 

Structures subjected to random excitations with uncertain 
system parameters degraded by surrounding environments (a 
random time history) are studied. Methods are developed to 
determine the statistics of dynamic responses, such as the 
time-varying mean, the standard deviation, the autocorrelation 
functions, and the joint probability density function of any 
response and its derivative. Moreover, the first-passage 
problems with deterministic and stationary/evolutionary 
random barriers are evaluated. The time-varying (joint) mean 
crossing rate and the probability density function of the first- 
passage time for various random barriers are derived. 

Introduction 

Many structures that require high performance or high 
reliability and durability often operate under complex 
environments including random excitations and random 
temperatures. These excitation and temperature variations 
not only degrade the material but also cause an additional 
randomness in the uncertain material behavior. In addition, 
a small variation in the structural shape may have an important 
effect on the structural responses. Therefore, the deterministic 


structural analysis is not sophisticated enough to quantify the 
uncertain structural responses, and the design requirements 
cannot be satisfied. For this reason a methodology has been 
developed for a probabilistic dynamic analysis of stochastic 
structures degraded by surrounding service environments. 
This methodology evolves from the concept of NESSUS 
(Numerical Evaluation of Stochastic Structures Under Stress), 
which is a probabilistic structural analysis code developed 
under NASA’s Probabilistic Structural Analysis Methods 
Project (Chamis, 1986). The methodology consists of five 
parts: (1) random process decomposition, (2) probabilistic 
material behavior model, (3) perturbed dynamic analysis of 
uncertain systems, (4) the first-order-second-moment method, 
and (5) a set of reliability algorithms. The random process 
decomposition procedure expresses the random phenomenon 
by a set of independent random variables and their respective 
characteristic functions. The probabilistic material behavior 
model defines the relationship between material properties 
(Young’s modulus, damping ratio, thermal expansion 
coefficient, strength, etc.) and the random phenomena 
(temperature, fatigue, etc). The perturbed dynamic analysis 
of uncertain systems produces perturbed time histories 
(dynamic responses), which are used to determine the response 
functions numerically in terms of independent random 
variables at each time. The first-order-second-moment 
method redefines the limit state function in the normalized 
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probabilistic space. A set of reliability algorithms arc then 
used to determine the important response statistics and to 
solve the first-passage problems. 

This newly developed method is demonstrated by performing 
a transient analysis of a stochastic structure subjected to random 
excitations under severe random temperature conditions. The 
structure is modeled by a single-degree-of-fireedom oscillator, 
or only the fundamental mode is considered in die dynamics 
analysis. The natural frequency of this oscillator, which is a 
function of the uncertain system parameters (geometry, material 
properties, etc.), is randomly degraded by the temperature 
process. The excitations can be either stationary or 
evolutionary random processes. Second-order statistics of 
any response, such as the time-varying mean, the standard 
deviation, and die autocorrelation functions, are determined. 
Also, the first-passage problems with constant or stationary 
evolutionary barriers are evaluated. A random barrier, such 
as the uncertain material strength or clearance for maximum 
displacement, is composed of its original randomness and of 
a variation due to the effect of random temperature and fatigue. 
The probability events in which crossing does or does not 
occur at any time are defined in terms of a reliability index 
and sensitivity factors from the first-order-second-moment 
method. The time-varying mean crossing rate and the joint 
mean crossing rate are derived. The probability density 
functions of the first-passage time for various random barriers 
are determined by two methods. The details of the method- 
ology are explained in this paper. The application of this 
methodology to a multiple-degree-of-freedom system in 
conjunction with NESSUS is also addressed. 


Random Process Decomposition 

For a random process in time, such as load and temperature, 
a Karhunen-Loeve expansion is used (Span os and Ghanem, 
1989). This expansion is a representation of the random 
process in terms of independent random variables U n and their 
associated time-varying characteristic functions f n (t), as shown 
in equation (1). 

S(0 = 5(0+i> n /„(0 0> 

n=l 


where S (t) is the mean value of the random process and r is 
the number of terms included. The smoother the random 
process is, the smaller the number r is. 


Probabilistic Material Behavior Model 

A generic material behavior model (multifactor interaction 
model) (Boyce and Chamis, 1988) is used to synthesize the 
variations of material properties for structures operating under 


hostile environments, such as high or cryogenic random 
temperature conditions. The fundamental assumption for this 
model is that the material behaviors can be synthesized by 
independent random variables V. The general form of this 
model is shown in equation (2), 


M P = 




( 2 ) 


where N is the number of environmental effects (temperature, 
fatigue, etc.); MpQ, Vp V 0 are independent random variables; 
Mpq is die reference material properties; VpVo, and Kare final, 
reference, and current environmental variables, respectively; 
and Mp is the current material property. The exponent a is 
determined from available experimental data or can be 
estimated from anticipated material behavior due to a particular 
environmental effect 


Perturbed Dynamic Analysis of 
Uncertain Systems 

The uncertain structural parameters, such as geometry or 
material properties, are random processes in space. The 
uncertain service environments, such as random excitations 
and temperature, are random processes in time. Furthermore, 
the material properties are affected and degraded by 
temperature. Because a proabilistic structural analysis by the 
first-order-second-moment method requires the input to be 
independent random variables, a decomposition procedure is 
applied to expand the random process in terms of independent 
random variables and their respective time-varying character- 
istic functions. Some other variables, such as initial boundary 
conditions, damping ratio, and random amplitude for harmonic 
loads, can also be included in the analysis. 

A dynamic analysis is performed repeatedly to generate a 
set of parallel time histories (dynamic responses) corresponding 
to a small change in each of the independent random variables. 
These perturbed time histories provide information on how 
the uncertainties in die problem parameters propagate through 
time and afreet the transient responses. The structural responses 
R(t) f in terms of independent random variables U it can be 
expressed as 


N N W 

RU) = a 0 + y £a,U i + 2 t b l U? 
l m 

where a j and b t are the coefficients computed by using 
perturbed time histories in a least-squares procedure. In the 
following, a first-order-second-moment method is used, 
together with the explicit response function determined by 
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equation (3), to derive the probabilistic information of 
structural responses. Also derived are the time-vaiying (joint) 
crossing rate and the probability density function of the first- 
passage time. 



Statistics of Dynamic Responses 

The dynamic responses of uncertain structures subjected to 
random excitation are also random processes in time. 
Although the first- and second-order statistical properties do 
not characterize a random process completely, they still 
contain the most important information of that process. 
Subsequent structural reliability assessments will be based 
on this information. In this paper the mean, the standard 
deviation, and the autocorrelation function of any random 
response process are derived in appendix A together with the 
joint probability density function (pdf) of any response and 
its derivative. The results are summarized in the following 
paragraphs. 

The cumulative distribution function (cdf) F, the mean 
E[R(t)] , and the standard deviation a for any structural response 
R(t) (displacement, velocity, acceleration, stresses, etc., at any 
time) are 

fy, )(*,> = *(-#) (4) 

_ k 

£[^(o]=/?(o=S 2 -^( z ‘) Az ( 5 ) 

i=l 

and 

( 6 ) 


where <X> is the standard normal cumulative distribution 
function; P is the reliability index from first-order-second- 
moment (FORM) analysis (Lind et al., 1985); k is the number 
of discretizations; and is the probability density function 
(pdf) of R(t). 

Defining X(t) as a random process (stress, displacement, 
etc.) and letting R 2 and R 2 be any two responses from X(/;), 

X(tj) t X(tj\X(t 2 ), X(t 2 ), Wf 2 ), the joint cumulative distribution 
function (cdf) F/^ derived in appendix A is 

1 - 1 ,...,* ;'= 1 / ( 7 ) 

where (i^yp are any pair of realizations; k and l are the number 
of discretizations for R j and R 2 , respectively; p,y is determined 
by equation (A 10); and 0 is the standard normal pdf. Letting 
R]~X(tj) and R 2 m X(t 2 ), the autocorrelation function of the 
random process X(t) is 



where fxqtfXQj is the joint pdf of X(tj) and X(t 2 ). 

Letting Rj m X(t) and R 2 m X (t), equation (7) becomes the joint 

cdf of X(t) and X (t). The mean rate at which a random process 
X(0 crosses a random barrier $(t) can be calculated by 

v*(0 = J~ f X(l )X( (9) 

where fx( t )x(t) “ d* J 0501 Pd* of *(# *(0- 

Probability Density Function of 
First-Passage Time 

The objective of a probabilistic dynamic analysis is to assess 
the reliability of structures subjected to random excitations. 
One of the failure modes concerning the structural reliability 
is so-called first-excursion failure. It states that the structure 
fails at the first excursion of a response to a given threshold. 
The response can be displacement, stress, strain, or eigenvalue. 
The threshold can be thought of as the maximum clearance 
for displacement at any time, the material strength to resist 
crack initiation caused by high stress, or the upper and lower 
bounds for eigenvalues within which responses can be 
minimized. The reliability for such problems is often 
represented by the probability density function of the time 
when the response fast crosses the threshold (pdf of the first- 
passage time). In this analysis the threshold is modeled as a 
deterministic time function and an evolutionary random 
process. The probability density function of the first-passage 
time is calculated by two methods, (1) a crossing-rate-based 
method and (2) an equivalent-system-based method, as 
described here. 

Crossing-Rate-Based Method 
The probability density function /j(r) of the first-passage 
time can be obtained by solving the integral equation 

/,(0+j**(f|T)/i(T)<ft = V*(0 (10) 

where £(f) is a random barrier process; vgt) is the mean 
crossing rate; Jf(/|T) is the crossing rate at t conditional on 
the first crossing having occurred at z. It is an integral 
identity expressing the fact that either a barrier crossing at 
time t must be the first or the first crossing must have occurred 
at some previous time z. Thus, the crossing rate can be given 
as the sum of two contributions from the left-hand side of 
equation (10). The simplest approximation of the pdf of the 
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first-passage time is from the approximation K(t\z) m vgt)\ that 
is, the crossing rate at t is assumed to be independent of a prior 
first crossing at t. In this case 

/ 1 (/)=v^(0- v 4(o£/i(t)A (ll) 


In the case where the Poisson hypothesis about independent 
barrier crossings is invalid, an improved approximation for 
the kernel AT((|t) is necessary. If the requirement that a barrier 
crossing at r should be first is dropped, then 


AT(r|x)= 


Vtt(*» T > 

v^(t) 


( 12 ) 


The pdf of the first-passage time is then calculated by die 
equation 


/.(0 = v 5 (r)-^-^j2/,(T)* (13) 


/l(«Af)Af = Prob 


= Prob 


= Ptob 


*i—i 


(ga^Oand ga_i>0)n(ft>0 and g;_i>0) 
1*1 


(gaSO)Dft >0 

r*l 


(MaSO)nV>0 

i*l 


(17) 


n-l 

The probability event Qa/, >0 is approximated by an 
i*l 

equivalent event A# e (it-1)>0, which is determined by a 
recursive equation as follows: 


A/ e (A-l)>0 = M ( *_ 1) >0fjM e (it-2)>0 Jk = 2 n 

(18) 


Therefore, 


where \)i*(r,T) is the joint mean crossing rate. Letting fn At 
and t * i At, equation (13) can be discretized as 

fl(nAt) = v^(nAf)- X ^ ( 14 > 

The mean crossing rate and the joint mean crossing rate are 
derived in appendixes B and C, respectively. 

Equivalent-System-Based Method 
The probability density function / j(n AO of the first-passage 
time is defined by the equation 

f\(n Af)A/ = Prob jpositive crossing in time[(/i - l)Ar,n A/] 
QNo crossing in time [(i-l)A/,i A/]i 

M J 

(15) 

Assume A/ is small enough to allow at most one crossing 
occurrence in this time interval. Defining the limit state 
function g i and its corresponding safety margin Af, by the 
equations 

N 

gi -^(iAt)-X(iAt) M i = fr-'£ 0 L ir U r i = l, .... n (16) 

r=l 

where £ is a random barrier process and X is a random response 
process; (3 and a are the reliability index and sensitivity factor, 
respectively; U r is a standardized normally distributed random 
variable; and N is the number of independent random variables. 
Substituting equation (16) into equation (15) gives 


fl(n At) At = Prob [ Af„ £ 0 and Af* (n - 1) > 0] (19) 

The probability cm the right-hand side of equation (19) will be 
calculated by equation (B8). 


Numerical Examples and Discussion 

The algorithms derived previously were applied to a single- 
degree-of-freedom oscillator shown in figure 1. The initial 
system parameter K 0 was affected by a stationary random 
temperature process 7X0, and die current system parameter 
K(t) was characterized by the following material behavior 
model: 


*<0 = *o 


Tf-TQ) 

t f -t 0 


( 20 ) 


where Kq Tp and T 0 are independent random variables and 
7(f) is a stationary random process whose covariance kernel is 
defined by the equation 

Cov [r(f),r(t)] = o* e" c *H (2i) 

where a T is the standard deviation of 7X0 and the constant Cj 
is the parameter that controls the smoothness of the random 
process. If Cf 0, the random process reduces to a single 
random variable. When Cj approaches infinity, the random 
process becomes a white noise. 

The natural frequency of the oscillator is obtained by the 
equation 
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where M is the lump mass. Substituting equation (20) into 
equation (22) gives 



Therefore, i o n (t) becomes a random process in time. The 
equation of the motion is then defined by the equation 

i'+2p© B (f)i+c)j;(r)jt = F(0 (24) 

where f) is the damping ratio. In the following numerical 
examples, these values are assumed: M • 300 000; 

ETKq] - 30 000 000; a^) - 600; E[T F ] « 2100; a(7>) - 42; 
T 0 -260 ; a - 0.25; £[7^)1 - 1800; afTO! - 36; C r « 0.25; 
P - 0.10. 

The oscillator is subjected to an evolutionary random 
excitation F(f), which is defined by 

F(t) = H(t)S(t) (25) 

where H(t) is a deterministic function defined by 

//(0=- t <t c 

= 1 t>t c (26) 

as shown in figure 2. Here 5(0 is a stationary random process 

with mean E[S(jt)] - 5(0 * 20 and standard deviation c s - 2. 
The covariance kernel of 5(0 is defined by the equation 

Cov[S(/),S(r)] = <7^' Cl H (27 ) 

Again, the constant C$ will determine the smoothness of 
the random excitation. 

A procedure is used next to decompose the random excitation 
5(0 &nd the random temperature 7\t) into a set of independent 
random variables and their respective time-varying character- 
istic functions. In the time interval [- / r /,] the process 5(0 can 
be expanded by the equation 

S(t) = ?(»)+ 5>/i( 0+ i>,/y(0 (28) 

i=0 j=Q 

where b i and bj are the zero-mean independent random 
variables and 



in which to,- and © y are the solution to the transcendental 
equations 

Cg - ©,• tan (©,/, ) = 0 (31) 

and 

(Oj + Cg tan (ci)y/<) = 0 (32) 

The random temperature 7(t) will be decomposed in a similar 
way. 

Statistics of Dynamic Responses 
The covariance of the random excitations is defined by 
equation (27) with C s - 5. This value will produce a 
nonsmooth random process with a correlation function similar 
to that for earthquake records. Thirty-six independent random 
variables are used to represent this process. The covariance 
of die random temperature is defined by equation (21) with 
C T m 0.25 because temperature oscillates slowly in time. Ten 
independent random variables are used to define this process. 
The mean, the standard deviation, and the autocorrelation 
function of the dynamic displacement were determined 
from equations (4) to (6). The results were compared with 
Monte Carlo simulation solutions and are shown in figure 3. 
Good agreement is observed. It is also noticed that in the 
time interval (0,2), the displacement process X(t ) is 
nonstationary with varying mean, standard deviation, and 
autocorrelation function. However, the process becomes 
stationary after 2 seconds. 

Probability Density Function of First-Passage Time 
The first-passage problems have been studied with different 
barriers and loading conditions. Two types of random barriers 
are used; they are defined by equations (33) and (34). 

Type A; 

*<0 = $ofl-7-l (33) 

V l F) 
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Type B: 



where £q is the reference barrier, t F is the final time, and 1 ftp 
is die degradation slope, which reflects the environmental 
degradation on the barriers. The type B random barrier includes 
the temperature effect more specifically. Two types of random 
loading are considered. The fust type is the (me with C 5 in 
equation (27) equal to 5. This type of loading is simil a r to 
earthquake records, which are nonsmooth random processes. 
The second type of loading has C 5 * 1, corresponding to a 
smooth random process. Two different degradation slopes 
are used, 0.02 and 0.15, to represent slow and fast 
environmental degradation. The mean reference barriers are 
chosen such that they correspond to o and 3o levels of the 
response process at t ■ 0. The statistics of random loads and 
barriers for each problem are listed in table I. 

The probability density functions of the first-passage time 
were calculated by both the crossing -rate-based method and 
the equivalent-system-based method. The results were 
compared with Monte Carlo simulation, as shown in figures 4 
to 7. With type A deterministic barriers, a nonsmooth random 
excitation (C s * 5), and a small barrier degradation slope 
( 0 . 02 ), the crossing-rate-based method performed well, as 
shown in figures 4(a) and (b), but the equivalent-system- 
based method did not. However, the equivalent-system -based 
method gave satisfactory results for problems with a large 
barrier degradation slope (0. 15), as shown in figures 4(c) and 
(d), but the crossing-rate-based method did not With the 
same barrier conditions as for problems in figures 4(a) and (b) 
and a smooth random excitation (C 5 * 1), the results from the 
equivalent-system-based method improved compared with 
figures 4 (a) and (b) as shown in figures 5(a) and (b). The pdf 
of the first-passage time predicted by the crossing-rate-based 
method for a large degradation slope also improved compared 
with figiffes 4(c) and (d), as shown in figures 5(c) and (d). 
Figure 6 shows the pdf of the first-passage time with nonsmooth 
random excitation (C 5 ■ 5) and type B random barriers. As 
expected, the equivalent-system-based method performed well 
for large degradation slopes, as shown in figures 6 (c) and (d). 
It also gave a satisfactory answer for a small degradation 
slope, as shown in figures 6 (a) and (b). For problems with 
type B random barriers and subjected to smooth random 
excitations (C 5 - 1 ), some improvements in the results 
predicted by the equivalent-system-based method were also 
observed for a small degradation slope ( 0 . 02 ), as shown in 
figures 7(a) and (b). As usual, it was good for a large 
degradation slope, as shown in figures 7(c) and (d). 

In general, the crossing-rate-based method is good for 
problems with a small degradation slope ( 0 . 02 ) subjected to 
nonwhite ( C s <* 5) excitations and for either low or high 
thresholds. The equivalent-system-based method is good for 
problems with a large degradation slope (0.15) subjected to 
nonwhite ( C s £ 5 ) excitations and for either low or high 


thresholds. However, when any of the following conditions 
are present, the results of the equivalent-system-based method 
will converge to the solutions by Monte Carlo simulation. 
The conditions are (1) smoother random excitations, ( 2 ) more 
uncertainties in the barrier, and (3) larger degradation slope. 
The reason is stated as follows: In order to have the equivalent- 
system-based method perform well, either the reliability index 
P; defined in equation (16) decreases rapidly or the correlations 
between the safety margins are large. Otherwise, a quadratic 
equivalent system should be developed to improve the 
accuracy. The degradation slope essentially controls the 
decreasing speed of The degree of uncertainty in the barriers 

and the smoothness of die random excitations, on the other 
hand, determine the correlations of the safety margins. 

Although the proposed method was applied to a single- 
degree-of-freedom oscillator, it can be directly applied to a 
multiple -degree -of-freedom system as addressed in the next 
section. The computing time required by Monte Carlo 
simulation is highly dependent on the degree of freedom of 
the finite element model. However, the computing time 
required for the proposed methods is dependent on the number 
of random variables and is not very sensitive to the degree of 
freedom of the finite element model. Therefore, they are 
more attractive for problems with complex structures or for 
nonlinear problems where analytical solutions do not exist It 
is necessary to point out that tire proposed method is not 
suitable for problems with white noise excitations because the 
infinite number of independent random variables required to 
represent this special random process makes the proposed 
method inefficient. 

Application to Multiple-Degree-of- 
Freedom System in Conjunction With 
NESSUS 

The algorithms derived in this paper were verified by 
applying them to a single-degree-of-freedom oscillator. 
However, the application to multiple-degree-of-freedom 
structures is straightforward and is discussed here. The 
equation of motion for such a system is defined by 

[M]x(/) + [c]x(r ) H-[A']x(f) = F{t ) (35) 

where the mass matrix [Af], the damping matrix [C], and the 
stiffness matrix [K] can be random processes in space. 
Because the explicit response functions in terms of 
independent random variables are required for the methods 
developed previously, a random process decomposition 
procedure is needed to represent the correlated random 
process in space by a set of independent random variables. 
This can be achieved by using NESSUS/PRE. NESSUS/ 
PRE is a preprocessor used for the preparation of the statistical 
data needed to perform the probabilistic structural analysis. 


6 


It allows the user to describe the uncertainties in the structural 
parameters (primitive random variables) at nodal points of a 
finite element mesh. The uncertainties in these parameters 
are specified over this mesh by defining the mean value and 
the standard deviation of the random variable at each point, 
together with an appropriate form of correlation. Correlated 
random variables are then decomposed into a set of indepen- 
dent vectors by a modal analysis. 

The perturbed time histories of any structural response 
corresponding to a small variation of each independent random 
variable can be generated by any finite element code with 
transient dynamic capability. The Finite Element Methods in 
NESSUS, called NESSUS/FEM, is one of the codes that can 
be used for this purpose. The explicit response functions can 
then be computed accordingly. 

NESSUS/FPI (Fast Probability Integrator) (Wu, 1985) 
provides the first-order-second-moment method necessary for 
determining the probabilistic responses and solving the first- 
passage problems. It is obvious that the newly developed 
methods not only utilize NESSUS but also enhance its 
capability to perfQrm various dynamic problems. It is worth 
pointing out that when the mass, damping, and stiffness 
matrices are independent of time, the transient dynamic 
responses can be obtained in closed forms by a modal 
analysis. The perturbed eigenvalues and corresponding 
eigenvectors are computed by using an efficient algorithm in 
NESSUS/FEM. 


Conclusions 


g‘=R{t)-i (Al) 

where R(t) is the structural response (X(r), X(r), X(t)) at 
time t, which is a function of independent random variables, 
and z is any realization. This limit state is subsequently 
transformed into a standardized normally distributed 
probability space and is replaced by a linear safety margin M 
through a fid-order-second-moment (FORM) analysis, where 
M is defined by 

W'=p'-£a' r f/ r (A2) 

r=l 

and where p* and a * are the reliability index and the sensitivity 
factor, respectively; U r is a standardized normally distributed 
random variable; and N is the number of independent random 
variables. The probability that 0 or R(t) £ z is determined 
by 

Probe's o)=Prob[/?(0< z] = Prob(M'£0) = <!>(-£') (A3) 

where O is the standard normal cumulative distribution 
function. The cumulative distribution function (cdf) F% for any 
structural response R(t) at any time can be generated from 
equation (A3) by selecting a number of realization z from an 
appropriate response range, as shown in equation (A4). 

/*(!)(*,) = *(-P') * = 1 * (A4) 


Methods have been developed for the probabilistic dynamic 
analysis of uncertain structures subjected to random excitations 
with random environmental effects on the uncertain structural 
parameters. The statistics of dynamic analysis such as time- 
varying mean, standard deviation, and autocorrelation function 
were determined. The first-passage problems, with barriers 
being deterministic functions and stationary/evolutionary 
random processes, were studied. The mean crossing rate and 
the joint mean crossing rate were derived. Two methods were 
proposed to determine the pdf of the first-passage time. The 
crossing-rate-based method is for problems of nonsmooth 
random excitations with small uncertainties in the barrier and 
small environmental effect on the barrier. The equivalent- 
system-based method is for problems of either large barrier 
degradation rate, smooth random excitations, or large 
uncertainties in the barrier. 


APPENDIX A 

STATISTICS OF DYNAMIC 
RESPONSES 


If fnt) I s denoted as the probabilty density function (pdf) of 
R(t), the mean and the standard deviation of R(t) can be 
calculated by the equations 

£[R(t)] = R(l) =X z </4)( 2 i) Az (AS) 

1=1 

and 

< y 4)=Ji[ z *-^)f/it( 2l )Az (A6) 


The joint cdf is derived as follows: Let Rj and R 2 be 

any two of X(t), X(x), X(t), X (x), X (/), X (x), etc., where X(t) 
is a random response process. Define the limit state functions 
S, and gj by 

gi= R l~ I i> gj = R 2~yj ( A7 > 


where z i and y. are different realizations. Also define their 
corresponding linear safety margins A/ ( and Mj by 

N N 

Mi = Pi riUr, Mj = P , -Xct rjV, (A8) 

r=l r= 1 


In order to calculate the response statistics, a limit state 
function g is first defined by 
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The joint cdf will be determined by the equation 

/ A#t»( z ^y) =Prob(/? i Sz * 81,(1 R 2-y^ 

= Prob(Afj 5 0 and Mj S 0) 

= 4>(-Pi) *(-P;) + (A9) 

1 = 1,..., k, j = 


Pij =X a n a /7 < A1 °) 

r= 1 

Let J7; - X(r) and R 2 m X(x) and calculate the autocorrelation 
function of X(t) by the equation 


££[*, ~ fjOjfe - 

CT X(r) CT X(-r) 

(All) 


where fx(t)X(i) is the joint pdf of X(t) and X(x). 


APPENDIX B 

MEAN CROSSING RATE v % (t) 


The corresponding linear safety margins M 1 and in the 
standardized normally distributed probability space are 

N N 

M 1 = p' - £a*C/ r , M ,+a * = P ,+Ar - £a‘ +A/ l/ r (B5) 

r=l r=l 


Substitute equations (B4) and (B5) into equation (Bl) to get 

= Prob(g‘>0 and 0) = Prob(Af'>0 and M ,+eu SO) 

(B6) 

Define the probability events A, B, A , and B by 

A=M l <,0 A=M'>0 

B=M , +*Z 0 F=Af' +A *>0 (B7) 

Then substitute equation (B7) into equation (B6) to get 

p;=Ptob(AB) 

= Prob(fl)-Prob(>4B) 

= Prob(M ,+ "S0)- ProbfAf SO and M' SO) 

= *(-?'“)-[*{-{*"“ )*(-/>' ) + £ j 

(B8) 


The mean rate at which a random process X(r) crosses a random 
barrier \{t) is determined by the equation 


^Average number of positive crossings in time (t,t + Ar)] 

= Prob[Positive slope crossings in time (f ,r + A/)] 

= Prob[x(r) < £(r) and X(r + Ar) > £(r + A/)] 

(Bl) 


'K 


where £(r) is the barrier. Therefore, 
v^(r)Ar = /’ c ' 


(B2) 


where 0 is the standard normal pdf and 

P =ix a r (B 9 ) 

r=l 


APPENDIX C 

JOINT MEAN CROSSING 

RATEv^(f,x) 

The joint mean rate at which a random process X crosses a 
random barrier \ at times t and X is derived by the following 
equations: 


Let / ■ n At, the crossing rate then becomes 

v t (»4/) = -£ (B3) 

where PL is determined as follows: The limit state functions 
g* and are defined by the equations 

g'=$(r)-X(r) g' +A, = ^ + A/)-X(/ + A/) (B4) 


[Average number of positive crossings in time(r,r + Ar) 
and positive crossings in time (t,t+A/)] 

= Probf Positive crossings in time (t,t + A/) and 
positive crossings in time (x, t + A/)] 

= Prob[X(r) < £(r) and X(r + Ar) > $(r +A/) and 
X(x)<£(x)and X(x+A/)>£(x+Ar)] 

= P" (Cl) 
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where t > T. Therefore, the joint mean crossing rate is 
determined by the following equation with t m n At and 
t - 1 An 

Vff(n &t,i Af)= (C2) 

ss (A/) 2 

where P£is determined by first defining the limit state functions 
gj, g} S3- “ d U 

gl =k(t)-X(t) g 2 =t,«+Al)-X(t+*t) 

«4 = ^(T+A/)-X(T+A/) (C3) 

The corresponding linear safety margins Mj, M 2 , M 3 , and 
In die standardized normally distributed probability space are 

N 

M i = fr-'£a ri U r i = 1,2, 3,4 (C4) 

r=l 

Then substituting equations (C3) and (C4) into equation (Cl) 
gives 

P”= Prob^g'>0 and g ,+A/ <oj and (g T > 0 and g 1+A, <ojJ 

= Prob^Af f >0 and Af ,+A, £o)and ^AfSOand M T+A, SojJ 

(C5) 

Also the probability events are defmed as 

A »Mj<0 ~A — A/j> 0 
B = M i < 0 B = M 2 > 0 

C = Af 3 <0 C = A/ 3 > 0 

D=M 4 < 0 D= M a > 0 (C6) 

Substituting equation (C6) into equation (CS) gives 

/>"= Prob( AS CD) 

= Prob(BD) + Prob (ABCD) - Prob(ABD) - Prob(£C£>) 

(C7) 

To alleviate the computational difficulty in equation (C7), 
two equivalent probability events E e j and £ f2 were derived to 
replace probability events AB and CD, respectively. Therefore, 
equation (C7) becomes 

P”= Prob(BD) + Prob(£ el E, 2 ) " Prob(£ el D) - Prob(fi£ e2 ) 

(C8) 


Each term cm the right-hand side of the equation can be easily 
determined from equation (A9). The equivalent systems 
and are determined as follows: Because the safety margins 
Afj and M 2 are highly correlated, the probability event AB can 
be replaced by an equivalent event (Goll witter and 

Rachwitz, 1983) in such a way that 

Prob[Afl] = Probes 0 and M 2 Z 0] 

= Probf^jS 0] = Prob[E e i] (C9) 

where M el is die equivalent linear safety margin defined by 

‘'Ur (CIO) 

r=l 


where 

P el = -<I>[Prob(A£)] (Cll) 


and 



The probability event CD is replaced by equivalent event E e2 
in a similar way. 
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Correlation function 




Figure 3. — Statistics of displacement in equation (24) (C$ in eq. (27), 5). 
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First-passage time, sec 


(c) Mean reference barrier, 0.37; (d) Mean reference barrier, 0.47; 

barrier degradation slope, 0.15. barrier degradation slope, 0.15. 


Figure 4. — Probability density function of first-passage time (C$ In eg. (27), 5; 
coefficient of variation of reference barrier, 0; type A deterministic barrier model). 
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— Equivaient-system-based method 

— Crossing -rate -based method 

Monte Carlo simulation 




(a) Mean reference barrier, 0.37; 
barrier degradation slope, 0.02. 


(b) Mean reference barrier, 0.47; 
barrier degradation slope, 0.02. 




First-passage time, sec 

(c) Mean reference barrier, 0.37; (d) Mean reference barrier, 0.47; 

barrier degradation slope, 0.15. barrier degradation slope, 0.15. 



Figure 5. — Probability density function of first-passage time (C$ In eq. (27), 1; 
coefficient of variation of reference barrier, 0; type A deterministic banler model). 
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— Equivalent-ays tem-based method 

— Crossing-rats -based method 
— — Monte Carlo simulation 



First -passage time, sec 


(c) Mean reference barrier, 0.50; (d) Mean reference barrier, 0.74; 

barrier degradation slope, 0.1 5. barrier degradation slope, 0.1 5. 

Figure 6.— Probability density function of first-passage time (Csin eq. 
(27), 5; coefficient of variation of reference barrier, 0.1; type B 
random barrier model). 
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Probability density 


Equivalent-systsm-besed method 

Crossing-rate-based method 

Monte Carlo simulation 




(c) Mean reference barrier, 0.58; (b) Mean reference barrier, 0.74; 

barrier degradation slope, 0.1 5. barrier degradation slope, 0.1 5. 

Figure 7. — Probability density function of first-passage time (C$in eq. (27), 1; 
coefficient of variation of reference barrier, 0.1; type B random barrier 
model). 
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